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Doped  barium  cerates  BaCeo.gLno.iC^  containing  earth-rare  dopants  with  different  ionic  radii,  Ln  =  La, 
Nd,  Sm,  Gd,  Yb,  Tb  and  Y,  have  been  investigated  as  candidate  materials  for  fuel  cells  and  other  electro¬ 
chemical  applications.  The  synthesis  of  these  materials  was  performed  using  a  precursor  method  based 
on  freeze-drying,  which  allows  a  precise  control  of  the  homogeneity  of  the  ceramic  powders.  Dense 
ceramic  pellets  were  obtained  at  1400  °C  under  identical  sintering  conditions.  The  microstructure  of  the 
ceramics  exhibits  similar  features  with  relative  density  higher  than  95%  and  the  grain  size  decreasing 
as  the  ionic  radius  of  the  dopant  decreases.  Impedance  spectroscopy  measurements  were  performed  to 
study  separately  the  different  contributions  to  the  total  conductivity.  The  bulk,  grain  boundary  and  total 
conductivities  depend  on  the  ionic  radius  of  the  dopant,  reaching  a  maximum  for  Gd-doped  samples  with 
a  value  of  0.02  S  cm-1  for  the  total  conductivity  at  600  °C. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Perovskite-type  oxides  based  on  BaCe03  exhibit  high  proton 
conductivity  at  intermediate  temperatures,  being  of  interest  for 
different  electrochemical  applications  [1-5].  Proton-conducting 
oxides  can  be  used  as  electrolyte  for  intermediate  temperature 
fuel  cells  (<600  °C)  [6-8],  with  several  advantages  over  their  oxide 
ion-conducting  counterparts.  The  oxide  ion-conducting  fuel  cell 
produces  water  vapor  on  the  anode  side,  which  lowers  the  cell 
voltage.  This  further  reduces  the  fuel  utilization  and  requires  con¬ 
siderable  fuel  circulation.  In  contrast,  a  proton-conducting  fuel  cell 
produces  water  at  the  cathode  side  while  the  fuel  at  the  anode 
stays  undiluted  by  water,  increasing  notably  the  fuel  utilization 
and  keeping  a  high  cell  voltage.  Furthermore,  BaCe03 -based  elec¬ 
trolytes  exhibit  higher  conductivity  than  commonly  used  oxide-ion 
electrolytes  based  on  Zr02,  Ce02  and  LaGa03  and  offer  the  addi¬ 
tional  possibility  of  reducing  the  operating  temperature  down  to 
500  °C  with  thin  electrolyte  layers  [9]. 

The  conducting  properties  of  the  acceptor-doped  BaCe03  have 
been  widely  investigated  [10-14].  Oxygen  vacancies  are  introduced 
into  perovskite  structure  by  partial  replacement  of  cerium  by:  Y3+, 
Yb3+,  Gd3+,  Nd3+,  etc.  [  1 5-21  ].  The  proton  conduction  in  these  mate¬ 
rials  occurs  either  via  the  hydration  of  oxygen  vacancies  after  the 
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material  is  exposed  to  water-containing  atmospheres  and/or  by 
reduction  of  oxide  ions  in  the  lattice  [11].  Proton  conductivity  is 
predominant  at  low  temperature,  whereas  oxide  ion  conductivity 
is  the  main  contribution  at  high  temperature  [19,20].  In  addition, 
n-type  and  p-type  electronic  conductivities  are  also  observed  at 
high  temperature,  limiting  the  application  of  cerates  as  electrolyte 
materials  above  700  °C. 

The  total  conductivity  of  BaCei_xMx03_5  depends  on  the  nature 
of  the  dopant  cation  M  and  its  content  x,  since  both  factors  alter  lat¬ 
tice  cell  parameters  and  degree  of  the  lattice  distortion,  which  affect 
the  mobility  of  ions.  Different  relationships  have  been  proposed  to 
describe  this  dependence  in  terms  of  a  critical  radius,  which  min¬ 
imizes  the  changes  in  the  lattice  parameters  of  the  pure  barium 
cerate  [15]. 

The  grain  interior  (bulk)  conductivity  of  ceramic  electrolytes 
often  determines  the  overall  electrical  properties  at  high  temper¬ 
atures;  however  the  grain  boundary  contribution  is  significant  in 
the  intermediate  and  low  temperature  range.  The  properties  of  the 
grain  boundary  are  known  to  be  strongly  affected  by  the  ceramic 
processing  conditions  and  microstructure  and  it  can  exert  a  pre¬ 
dominant  effect  on  the  properties  of  the  ceramic  materials.  Thus, 
the  optimisation  of  ceramic  materials  should  not  only  be  based 
on  the  search  of  the  most  promising  phase  composition,  but  also 
requires  a  detailed  study  of  the  relationship  between  ceramic  pro¬ 
cessing  conditions  and  its  microstructure. 

The  main  synthesis  route  used  in  the  preparation  of  BaCe03  is 
the  solid  state  ceramic  method,  which  requires  firing  temperatures 
as  high  as  1600  °C  to  obtain  both  a  single  phase  and  dense  samples 
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Table  1 

Structural  and  microstructural  parameters  for  BaCeo.9Gdo.i03_5  series. 


Ln 

v(A3) 

Ri  (A) 

Tolerance  factor  (t) 

Free  volume  Vf  (A3) 

Dg(^m) 

Relative  density  (%) 

La 

341.90 

1.032 

0.9307 

30.59 

1.5 

95 

Nd 

341.25 

0.983 

0.9329 

30.45 

1.4 

96 

Sm 

341.15 

0.958 

0.9340 

30.44 

1.3 

100 

Gd 

340.80 

0.938 

0.9348 

30.36 

1.2 

98 

Y 

340.23 

0.900 

0.9363 

30.23 

1.1 

95 

Yb 

340.80 

0.868 

0.9377 

30.21 

1.0 

98 

Tb4+,3+ 

339.40 

- 

- 

- 

2.0 

91 

[21].  This  method  has  several  disadvantages,  such  as  formation  of 
undesirable  phases,  poor  chemical  homogeneity  and  large  grain 
ceramic  sizes,  which  negatively  affect  the  mechanical  performance. 
In  addition,  evaporation  of  BaO  has  also  been  detected  when  the 
sintering  temperature  exceeds  1500°C  [22,23],  thus  reducing  the 
sintering  temperature  is  particularly  important  for  these  materi¬ 
als. 

An  alternative  preparation  method  based  on  freeze-drying 
precursor  route  was  recently  used  to  obtain  powders  of 
BaCeo.gYo.i  03_5  [23  ].  This  method  has  proved  to  be  very  efficient  for 
obtaining  stoichiometrically  controlled  polymetallic  compounds 
[24,25],  yielding  homogeneous  and  very  fine  powder  that  con¬ 
tributes  to  avoid  heterogeneous  distribution  of  cations.  Moreover, 
the  differences  in  transport  properties  related  to  sample  prepa¬ 
ration  can  be  avoided  by  using  identical  powder  preparation  and 
densification  processes. 

In  the  present  study,  the  perovskite-type  oxide  solid  solu¬ 
tions  BaCei_xLnx03_5  with  trivalent  dopants  of  different  ionic  size 
(Ln  =  La,  Nd,  Sm,  Gd,  Yb,  Tb  and  Y)  and  with  fixed  amount  of  their 
content  (x  =  0.1)  were  synthesized.  The  microstructure,  bulk  and 
grain  boundary  contributions  to  the  total  conductivity  were  stud¬ 
ied  under  different  atmospheres  with  the  emphasis  on  a  possible 
relationship  between  processing  and  microstructure. 

2.  Experimental 

2.1.  Synthesis 

Powders  of  BaCe!_xLnx03_5  (Ln  =  La,  Nd,  Sm,  Gd,  Y,  Yb  and  Tb) 
were  prepared  by  a  freeze-drying  method.  The  reactants  were 
high  purity  nitrates:  Ba(N03)2  (+99%  Aldrich)  and  Ln(N03)3-6H20 
(99.99%  Aldrich).  These  nitrates  were  previously  studied  by 
thermogravimetric  analysis  to  determine  the  correct  cation  com¬ 
position.  Stoichiometric  quantities  of  reactants  were  dissolved  in 
distilled  water  and  ethylenediaminetetraacetic  acid  (EDTA)  (99.7% 
Aldrich)  was  added  in  a  1:1  molar  ratio  ligand:metal  as  complex 
agent,  which  is  necessary  to  prevent  precipitation.  A  typical  prepa¬ 
ration  of  100  ml  contains  15  mmol  of  Ba2+.  The  pH  was  adjusted  at 
7  by  adding  ammonia.  The  solutions  were  flash-frozen  by  drop¬ 
ping  them  into  liquid  nitrogen  in  order  to  maintain  the  cation 
homogeneity  of  the  starting  solution,  and  then  dehydrated  using 
a  HetoLyolab  freeze-dryer  for  3  days.  In  this  way  dried  and  amor¬ 
phous  precursor  powders  were  obtained.  These  precursors  are  very 
hygroscopic  so  that  they  were  immediately  calcined  at  300  °C  to 
prevent  rehydration  and  to  produce  the  combustion  of  organic 
material.  The  precursor  powders  obtained  in  this  way  were  cal¬ 
cined  between  300  and  1000°C  for  1  h,  in  order  to  study  the  phase 
formation  behaviour  with  temperature.  All  the  thermal  treatments 
were  carried  out  in  alumina  crucibles  to  prevent  reaction  of  barium 
with  platinum  crucibles. 

The  as-prepared  oxide  powders  at  1000°C,  were  pressed 
at  200  MPa  into  disks  of  10  mm  diameter  and  1.5  mm  thick¬ 
ness.  All  the  pellets  were  sintered  in  air  at  1400  °C  for  4h. 
The  relative  density  of  the  sintered  pellets  was  calculated 
from  the  mass,  volume  of  the  sample  and  theoretical  density 


obtained  from  the  structural  analysis  by  Rietveld  refinement 
(Table  1). 

2.2.  Structural  and  microstructural  characterisation 

Room  temperature  X-ray  diffraction  patterns  (XRD)  were  col¬ 
lected  with  a  PANalytical  X’Pert  diffractometer  with  Cu  Ka 
radiation  and  the  X’Celerator  detector.  The  scans  were  performed 
in  the  20  range  (15-100°)  with  0.016°  step  for  2h.  Structure 
refinements  were  performed  using  the  FullProf  [26]  and  the  ICSD 
database  for  the  structural  models  [27]. 

X-ray  photoelectron  spectroscopy  (XPS)  spectra  were  recorded 
on  a  Physical  Electronics  PHI-5700  spectrometer  equipped  with  Mg 
Ka  and  Al  Ka  X-ray  sources  and  a  hemispherical  electron  analyser. 
The  data  were  analysed  using  XPSpeak  4.1  software  [28].  The  exper¬ 
imental  peaks  were  deconvoluted  into  components  using  mixed 
Gaussian-Lorentzian  functions. 

The  microstructure  of  the  sintered  pellets  was  observed  by  Scan¬ 
ning  Electron  Microscope  (SEM)  using  a  Jeol  Ltd.,  JSM-6300.  All 
samples  were  covered  with  a  sputtered  thin  film  of  gold  to  avoid 
charging  problems  and  to  obtain  better  image  definition.  Aver¬ 
age  grain  size  of  sintered  pellets  was  estimated  from  the  SEM 
micrographs,  using  the  linear  intercept  method  from  at  least  50 
randomly  lines  and  two  different  micrographs  with  the  help  of  an 
image-analysis  software  [29]. 

2.3.  Electrical  characterisation 

The  electrical  characterisation  of  the  ceramic  material  was  car¬ 
ried  out  by  impedance  spectroscopy  under  different  dry /wet  gases 
(air  and  5%  H2-Ar).  Gases  were  humidified  using  a  gas-washer  at  a 
temperature  of  20  °C  to  ensure  constant  water  content  of  about  2%. 
Pt  electrodes  were  deposited  on  each  side  of  the  sintered  pellet  with 
Pt-ink  (Metalor)  and  then  thermally  treated  at  800  °C  for  30  min 
to  ensure  maximum  conductivity  and  adherence.  Impedance  spec¬ 
tra  were  obtained  using  a  frequency  response  analyser  (Solartron 
1260)  in  the  0.1  Hz  to  1  MHz  frequency  range  with  an  ac  signal 
of  50-1 00  mV.  The  spectra  were  acquired  during  the  cooling  pro¬ 
cess  between  800  and  1 50  °C  with  a  dwell  time  of  60  min  between 
consecutive  measurements.  Equivalent  circuits  were  used  to  fit  the 
spectra  and  to  study  separately  the  different  resistive  contributions 
(bulk  and  grain  boundary)  to  the  total  conductivity. 

Conductivity  dependence  as  a  function  of  oxygen  partial  pres¬ 
sure  (p02)  from  0.21  to  ~10-20atm  was  obtained  with  an 
electrochemical  cell,  which  consist  of  a  YSZ  tube  closed  at  one  end, 
with  an  electrochemical  pump  and  an  oxygen  sensor.  Before  tak¬ 
ing  measurements,  the  pellets  were  kept  at  700-800  °C  for  5  h  while 
the  system  was  flushed  with  humidified  5%  H2-Ar  gas  mixture.  The 
oxygen  partial  pressure  was  controlled  using  the  electrochemical 
pump  with  a  dc  source  (Yokogawa  7651)  and  the  oxygen  partial 
pressure  inside  the  cell  was  continuously  monitored  by  means 
of  the  oxygen  sensor  (voltage  read  by  Keithley  2700  multime¬ 
ter)  placed  close  to  the  sample.  The  samples  were  equilibrated  at 
each  oxygen  partial  pressure  at  least  for  2  h  before  acquiring  the 
impedance  spectra. 
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Fig.  1.  XRD  powder  patterns  of  BaCeo.gLno.i03_5  (Ln  =  La,  Nd,  Sm,  Gd,  Tb  and  Yb) 
prepared  by  freeze-drying  method  at  1000  °C  for  1  h. 


3.  Results  and  discussions 

3.1.  Structure 


The  perovskite  phase  crystallizes  after  one  single  thermal  treat¬ 
ment  of  the  precursor  powders  at  800  °C,  although  a  small  amount 
of  BaC03  impurity  is  still  detected  by  means  of  high  resolution  XRD. 
This  impurity  disappears  after  calcination  at  1000°C  for  1  h.  The 
XRD  patterns  corresponding  to  as-prepared  powders  at  1000°C 
(Fig.  1)  indicate  that,  in  all  samples,  the  perovskite  structure  is 
obtained  and  no  additional  phases  are  observed. 

In  order  to  improve  crystallization  of  the  compounds  and  facil¬ 
itate  their  structural  characterisation,  the  polycrystalline  powders 
were  pressed  into  pellets,  calcined  at  1400  °C  for  4  h  and  then  slowly 
cooled  in  air  to  room  temperature.  The  XRD  patterns  were  anal¬ 
ysed  by  the  Rietveld  method  using  the  structural  models  previously 
proposed  in  the  literature  (ICSD  72768)  in  the  orthorhombic  Pmcn 
space  group. 

The  lattice  cell  volume  was  plotted  against  the  ionic  radii  [30] 
of  the  six  coordinated  rare-earth  dopant  cation  (Fig.  2a).  In  this  fig¬ 
ure  a  linear  dependence  is  observed,  which  is  the  typical  behaviour 
expected  for  a  solid  solution.  Flowever,  this  trend  is  not  followed 
by  Tb-doped  sample,  if  the  perovskite  structure  is  supposed  to 
contain  only  six  coordinated  Tb3+  with  ionic  radius  of  0.90  A.  This 
sample  was  thermally  treated  under  a  continuous  flux  of  wet 
diluted  hydrogen  (5%  FI2-Ar)  at  800  °C  and  then  analysed  by  XRD. 
The  cell  volume  changed  significantly  from  337.56  to  339.40  A3 
after  the  thermal  treatment  in  oxidising  and  reducing  atmospheres 
respectively  (Fig.  2a);  however,  this  important  volume  change  was 
not  observed  in  the  other  samples.  These  results  seem  to  indi¬ 
cate  that  Tb4+  and  Tb3+  ions  coexist  simultaneously  in  reduced 
BaCe0.9Tb0.iO3_5  sample.  This  issue  was  further  confirmed  by  XPS 
and  impedance  spectroscopy  analysis  in  the  next  sections. 

It  is  generally  assumed  that  lower  distortion  of  the  lattice  cell 
and  larger  free  volume  favour  the  ionic  conduction  in  perovskites. 
In  order  to  correlate  the  structure  with  the  transport  properties,  the 
tolerance  factor  and  crystallographic  free  volume  were  determined. 


Fig.  2.  (a)  Variation  of  the  unit  cell  volume  and  (b)  free  volume  ( Vf)  and  the  tolerance 
factor  (t)  in  BaCeo.9Ln0.i03_3  (Ln  =  La,  Nd,  Sm,  Gd,  Tb,  Yb  and  Y)  compounds  as  a 
function  of  the  rare-earth  ionic  radius.  Open  circles  in  (a)  represent  the  cell  volume 
for  Tb-doped  sample  considering  only  Tb4+  or  Tb3+  in  the  structure. 


The  degree  of  distortion  of  AB03  perovskite  structure  is  given 
by  the  Goldschmidt  tolerance  factor  (t),  defined  as  [31  ]: 


rA  +  rO 
V2(Fb  +  7q) 


(1) 


where  rA  and  rB  are  the  ionic  radii  of  the  cations  occupying  the 
A  and  B  site  respectively,  while  r0  is  the  oxygen  ionic  radii.  The 
tolerance  factor  of  undoped  barium  cerate  is  0.9376  and  it  decreases 
linearly  with  increasing  the  earth-rare  ionic  radius  of  the  dopant 
from  0.9377  for  Yb3+  to  0.9307  for  La3+  (Fig.  2b).  Thus,  the  lattice 
distortion  of  BaCe03  perovskite  decreases  as  the  ionic  radius  of  the 
dopant  in  the  B-position  is  getting  smaller. 

On  the  other  hand,  the  free  volume  (Vf)  is  defined  as  the  lattice 
unit  cell  volume  minus  the  ionic  volume  of  the  constituent  ions.  As 
can  be  observed  in  Table  1  and  Fig.  2b  the  free  volume  decreases 
with  decreasing  the  ionic  radius  of  the  dopant  from  30.59  A3  for 
La3+  to  30.21  A3  for  Yb3+. 

Thus,  dopants  with  smaller  ionic  radii  exhibit  higher  toler¬ 
ance  factor  (increasing  the  ionic  mobility)  and  lower  free  volume 
(decreasing  ionic  mobility);  and  a  contrary  trend  is  observed  for 
dopant  with  larger  ionic  radii.  Flence,  cations  with  intermediate 
values  of  free  volume  and  tolerance  factor  are  expected  to  be  the 
best  ones  for  optimizing  the  ionic  conductivity  in  barium  cerates. 


3.2.  XPS  analysis 


Several  sintered  pellets  of  BaCe0.gLn0.i  03_5  were  exposed  to  oxi¬ 
dising  and  reducing  atmospheres  at  800  °C  under  a  constant  flow  of 
wet  air  and  wet  5%  FI2-Ar  gas  mixture  for  12  h.  The  pellet  surfaces 
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were  ground  with  SiC  paper  before  XPS  measurements  to  remove 
the  superficial  layer,  because  BaCe03  materials  are  susceptible  to 
surface  carbonation  in  the  presence  of  C02. 

It  should  be  commented  that  the  spectra  of  reduced  Gd  and  Yb- 
doped  samples  are  shifted  to  higher  binding  energies  compared 
to  the  oxidised  samples;  and  a  reverse  tendency  is  observed  for 
Tb-doped  samples  (Fig.  SI  and  S2,  Supplementary  Materials).  This 
behaviour  may  be  explained  by  the  different  conductivity  values 
of  these  samples  in  reducing  and  oxidising  atmospheres  as  fur¬ 
ther  discussed  below.  Gd-  and  Yb-doped  samples  exhibit  larger 
conductivity  values  under  oxidising  conditions  due  to  the  pres¬ 
ence  of  a  significant  p-type  contribution.  Conversely,  Tb-doped 
samples  have  very  low  conductivity  under  oxidising  condition.  In 
this  sample,  Tb  is  in  an  oxidation  state  +4  in  air  and  it  partially 
reduces  to  Tb3+  under  hydrogen  atmosphere,  increasing  the  oxy¬ 
gen  vacancy  concentration  and  enhancing  considerably  the  ionic 
conductivity. 

The  Ba  3d  core  level  has  a  quite  simple  spectrum  in  all  the  sam¬ 
ples,  a  typical  3d  doublet  with  FWHM  of  2.0-2. 1  eV  and  only  one 
contribution  (Fig.  SI  Supplementary  Materials).  The  Ba  3d5/2  peak 
shows  a  maximum  at  ~785  eV,  which  is  higher  than  expected.  This 
is  probably  related  to  the  low  electronic  conductivity  of  these  sam¬ 
ples.  Thus,  a  shift  in  the  binding  energies  of  ~4eV  was  applied 
in  order  to  ascribe  the  different  contributions  of  the  core  levels 
spectra. 

C 1  s  spectra  show  two  separated  peaks  with  at  least  five  different 
contributions  in  all  the  samples  (Fig.  S2  Supplementary  Materials). 
They  have  been  ascribed  to  adventitious  carbon,  different  C-0  and 
C=0  species  and  carbonates,  being  those  corresponding  to  C-0  and 
carbonates  the  most  intense  ones. 

The  spectra  corresponding  to  the  4d  core-level  of  the  rare-earth 
dopants  exhibit  very  low  intensity.  The  characteristic  4d  doublet 
of  the  rare-earth  dopants  is  difficult  to  distinguish  due  to  the 
low  fraction  of  these  cations  in  the  samples  and  their  very  low 
intensity  in  the  XPS  spectra.  A  comparison  between  the  spectra 
of  reduced  and  oxidized  samples  does  not  reveal  appreciable  dif¬ 
ferences  with  the  exception  of  the  Tb  4d  core  level.  In  this  case, 
the  spectrum  of  the  reduced  sample  shows  a  tiny  shoulder  at  low 
energy  side  compared  to  the  oxidised  sample,  which  is  associated 
to  the  presence  of  mixed  oxidations  states  of  Tb4+  and  Tb3+  ions 
(not  shown). 

O  Is  core  levels  show  a  rather  complicated  structure,  quite 
similar  in  all  the  samples.  Four  different  contributions  have  been 
considered  to  adequately  fit  the  experimental  spectra  (Fig.  3).  The 
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Fig.  3.  XPS  spectra  for  the  0  Is  core  level  of  BaCe0.9Gd0.iO3_5. 
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Fig.  4.  XPS  spectra  for  the  Ce  3d  core  levels  of  BaCeo.gGdo.iO^  (BCG)  and 
BaCeo.gTbo.i  03_,$  (BCT)  samples  as-prepared  in  air  and  after  thermal  treatment  under 
5%  H2-Ar  flow  at  800  °C. 

lowest  binding  energy  contribution  at  529.9  eV  might  be  assigned 
to  different  superficial  Ce-containing  oxides  and  hydroxides  [32]. 
The  second  contribution  is  associated  to  superficial  oxygen  contam¬ 
ination,  such  as  specimens  containing  C-0  bonds.  The  most  intense 
contribution,  about  45%  of  the  total  peak  intensity,  with  a  binding 
energy  higher  than  those  corresponding  to  simple  metal  oxides 
is  ascribed  to  the  oxygen  belonging  to  the  perovskite  structure. 
Finally,  a  small  contribution  with  less  than  10%  of  the  total  peak 
intensity  is  observed  at  higher  binding  energies,  which  is  probably 
related  to  Ln-0  bonds  of  the  trivalent  dopant  (Fig.  3). 

The  Ce  3d  core  level  data  show  the  usual  multicomponent  spec¬ 
tra  due  to  the  presence  of  different  Ce  oxidation  states  (Fig.  4). 
One  would  expect  negligible  Ce3+  contribution  in  oxidising  atmo¬ 
spheres,  however  partial  reduction  of  cerium  during  recording  of 
the  XPS  spectra  at  the  surface  of  the  pellets  is  possible  as  has  been 
previously  observed  [33  ].  It  is  worth  mentioning  that  all  the  spectra 
look  very  similar,  except  that  corresponding  to  the  oxidised  Tb- 
doped  sample  (Fig.  4),  which  exhibits  the  highest  content  of  Ce4+. 
This  seems  to  indicate  that  oxidised  BaCeo.gTbo.i  03_5  without  oxy¬ 
gen  vacancies  in  the  structure  shows  lower  degree  of  reducibility 
than  those  containing  oxygen  vacancies.  It  should  be  noted  that 
a  similar  behaviour  has  been  also  observed  in  undoped  and  doped 
Ce02  [34].  Considering  also  that  Ce  cations  have  different  symmetry 
environments  due  to  the  presence  of  several  Ce-containing  com¬ 
pounds  at  the  surface  of  the  samples,  such  as  hydroxides,  there 
will  be  more  than  one  contribution  to  the  spectra  for  each  oxida¬ 
tion  state.  Thus,  the  Ce3+/Ce4+  fraction  is  difficult  to  determine  in 
these  compounds. 

3.3.  Ceramic  micro  structure 

The  microstructure  of  the  sintered  pellets  was  analysed  as  a 
function  of  the  dopant  type  from  SEM  micrographs.  SEM  images 
for  the  different  pellets  sintered  at  1400  °C  for  4h  are  shown  in 
Fig.  5.  No  appreciable  traces  of  impurities  or  phase  segregations  at 
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Fig.  5.  SEM  images  of  BaCeo.gLn0.i 03_,$,  Ln  =  (a)  La,  (b)  Nd,  (c)  Sm,  (d)  Gd,  (e)  Tb  and  (f)  Yb,  pellets  sintered  at  1400 °C  for  4  h. 


the  grain  boundaries  were  found  in  any  of  the  samples  studied.  The 
microstructure  is  well  developed  with  well  connected  grains  and 
relative  density  higher  than  95%  (Table  1).  Therefore,  the  possible 
effects  of  the  porosity  on  the  bulk  and  grain  boundary  conduction 
are  ruled  out.  In  general,  the  different  dopants  do  not  affect  signifi¬ 
cantly  the  microstructure  of  the  pellets  and  differences  are  mainly 
observed  in  the  average  grain  size.  The  undoped  BaCe03  sample 
exhibits  an  average  grain  size  of  2.5  |jim  after  sintering  at  1400  °C 
for  4h  [23],  however  smaller  grain  sizes  are  observed  in  the  doped 
samples,  which  is  usually  ascribed  to  lower  diffusion  across  the 
grain  boundary  during  the  sintering  process.  The  densification  pro¬ 
cess  occurs  in  a  first  step  and  then  further  grain  growth  is  observed. 
Both  processes  are  generally  controlled  by  grain  boundary  diffusion 
rather  than  bulk  diffusion.  An  enrichment  of  the  trivalent  dopant 
at  the  grain  boundary  region,  according  to  the  space  charge  model 
[35],  generates  a  concentration  gradient  between  the  grain  interior 
and  grain  boundary  that  decreases  the  mobility  of  cations  along  the 
grain  boundary  and  therefore  the  densification  and  grain  growth 
rate  decreases.  A  similar  effect  has  previously  been  observed  in 
other  materials,  such  as  doped-Ce02  [36]. 


The  variation  of  the  average  grain  size  with  the  ionic  radii  of 
the  rare-earth  dopant  shows  a  linear  behaviour  (Fig.  6),  which  is 
almost  identical  to  that  observed  in  the  lattice  cell  volume  (Fig.  2a). 
The  ceramic  grain  size  becomes  smaller  as  the  ionic  radii  of  the 
dopant  decreases  from  1.5  [xm  for  La3+  to  1  [xm  for  Yb3+  (Table  1). 
Flowever,  this  linear  relationship  does  not  hold  for  Tb-doped  sam¬ 
ple  sintered  in  air,  which  exhibits  the  largest  grain  size  ~2  fxm  and 
comparable  to  the  undoped  one.  This  effect  can  be  explained  taking 
into  account  that  Tb-doped  samples  sintered  in  air  do  not  contain 
oxygen  vacancies  in  their  structure,  because  Tb  is  in  oxidation  state 
+4,  and  therefore  the  densification  process  seems  to  be  similar  to 
that  occurring  in  undoped  BaCe03. 

3.4.  Electrical  characterisation 

The  impedance  spectra  of  BaCe0.gLn0.iO3_5  recorded  in  air 
show  similar  features  (Fig.  7)  with  three  separated  contributions 
ascribed  to  grain  interior  (bulk),  grain  boundary  and  electrode 
processes.  Similar  spectra  were  obtained  in  wet  5%  H2-Ar  gas 
mixture.  The  results  suggest  an  ideal  equivalent  circuit  consisting 
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Fig.  6.  Variation  of  the  average  grain  size  of  the  sintered  pellets  of  BaCe0.9Ln0.iO3_5 
(Ln  =  La,  Nd,  Sm,  Gd,  Tb,  Yb)  series  as  a  function  of  the  rare-earth  ionic  radius. 


of  (RC)  elements  in  series,  although  the  spectra  show  somewhat 
depressed  arcs.  Therefore,  the  alternative  equivalent  circuit  used 
was  (RbQbXftgbQgbXQe)  where  the  capacitors  C  were  replaced  by 
constant  phase  elements  Qwith  impedance  value  ofZ(±=  l/[QC/w)n] 
to  achieve  better  results.  The  following  parameters  were  fitted  for 
each  depressed  arc:  the  resistance  R2,  the  pseudocapacitance  Qi  and 
the  exponential  parameter  nz  of  bulk,  grain  boundary  and  electrode 
processes.  These  parameters  are  related  to  the  angular  relaxation 
frequency  w  and  true  capacitance  C  by  the  following  relation  [37]: 

wf  =  (RjC,)-1  =  (2) 

Bulk  and  grain  boundary  contributions  show  typical  capacitance 
values  of  pF  cm-1  and  nF  cm-1 ,  respectively.  Tb-doped  sample  dis¬ 
plays  different  impedance  spectra  in  air  (Fig.  7f)  with  only  one  broad 
semicircle  and  a  capacitance  value  of  6  pF  cm-1 ,  which  is  associated 
to  bulk  contribution.  Nevertheless,  the  spectra  of  this  sample  in  H2- 
containing  atmosphere  are  very  similar  to  those  of  the  rest  of  the 
samples  with  three  well  separated  contributions  and  much  lower 
resistance  due  to  Tb4+  ->  Tb3+  reduction  and  the  consequent  forma¬ 
tion  of  oxygen  vacancies  that  enhances  the  ionic  conductivity. 

Grain  boundary  and  grain  interior  processes  were  studied 
separately  in  the  low  temperature  range  (T<250°C)  were  both 
processes  are  discernible.  Conversely,  electrode  processes  are  dom¬ 


inant  at  higher  temperatures  (T>  400  °C)  and  only  total  conductivity 
was  obtained.  In  this  case,  the  equivalent  circuit  used  to  fit  the 
spectra  is  denoted  as:  FRs(^eiQeiX^e2Qe2)*  where  L  is  an  autoin¬ 
ductance  element,  which  is  associated  to  the  autoinductive  effects 
introduced  by  the  equipment,  and  Rs  is  a  serial  resistance  equal 
to  total  electrolyte  resistance;  Rei  and  Qei  are  the  resistance  and 
pseudocapacitance  of  each  electrode  process. 

The  relative  permittivity  of  the  bulk  was  estimated  from  the 
capacitance  values  by  using  the  following  relation: 


(3) 


where  L  and  S  are  the  thickness  and  electrode  area  of  the  pellets 
respectively  and  s0  is  the  vacuum  permittivity.  A  relative  permit¬ 
tivity  sb  ~  40  was  obtained,  which  is  comparable  to  the  previously 
reported  one  £b  ~  37  [38]. 

The  grain  interior  Rb  and  grain  boundary  Rgb  resistances  were 
used  to  obtain  the  bulk  and  macroscopic  grain  boundary  conduc¬ 
tivity,  taking  into  account  the  sample  geometry: 


(4) 


The  values  of  conductivity  were  plotted  using  the  Arrhenius 
equation: 


The  temperature  dependence  of  grain  interior  conductivity  for 
the  different  samples  in  wet  air  and  wet-5%  H2-Ar  gases  are  shown 
in  Fig.  8a  and  b  respectively.  As  can  be  observed  the  bulk  con¬ 
ductivity  depends  on  both  the  ionic  radii  of  the  dopant  and  the 
atmosphere  with  slightly  larger  values  of  conductivity  in  wet-5% 
FI2-Ar.  The  activation  energy  for  bulk  conduction,  in  the  tem¬ 
perature  range  of  1 50-250  °C,  decreases  slightly  with  decreasing 
ionic  radii  of  the  dopant  from  0.53  eV  for  La3+  to  0.46  eV  for  Tb4,3+ 
(Table  2).  In  addition,  the  dependence  of  the  bulk  conductivity 
as  a  function  of  the  rare-earth  dopant  size  shows  a  maximum 
for  Gd-doped  samples  (Fig.  9).  As  aforementioned,  the  best  bulk 
ionic  conduction  in  doped  perovskite  electrolytes  corresponds  to 
minimum  lattice  distortion  and  largest  lattice  free  volume.  The 
lattice  distortion  in  doped  BaCe03,  defined  by  the  tolerance  fac¬ 
tor,  is  lower  as  the  ionic  radius  of  the  dopant  decreases,  favouring 
the  ionic  conduction.  Conversely,  the  lattice  free  volume  decreases 
with  decreasing  ionic  radii  of  the  dopant  and  thus  reducing  the  ionic 
mobility.  Both  phenomena  are  involved  in  the  bulk  conduction  and 
they  can  explain  the  maximum  of  the  bulk  conductivity  observed 


Fig.  7.  Impedance  spectra  for  BaCeo.9Ln0.i03_5  pellets  sintered  at  1400  °C  for  4  h:  Ln  =  (a)  La,  (b)  Nd,  (c)  Sm,  (d)  Gd  and  (e)  Yb  measured  at  125  °C  and  f)  Tb  at  600  °C.  The  solid 
line  is  curve  obtained  by  fitting  the  data  with  equivalent  circuits. 
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Fig.  8.  Arrhenius  plots  of  the  grain  interior  (bulk)  conductivity  for  BaCeo.9Ln0.i03_5  samples  under  (a)  wet  air  and  (b)  wet  5%  H2-Ar  gases. 


Table  2 

Conductivity  and  activation  energies  for  the  bulk  and  grain  boundary  contributions  obtained  from  the  impedance  spectra  in  air  for  BaCeo.gLno.iO^  series. 


Ln 

0bUik  (fiScrrr1)200"c(air 

<Tgb  (|JiS  cm4 ),  200°c/air 

Fbuik  (eV)air 

Fgb  (eV)air 

<5gb  (nm)air 

La 

33 

0.04 

0.53 

0.82 

4.5 

Nd 

36 

0.22 

0.53 

0.75 

4.6 

Sm 

151 

0.41 

0.49 

0.76 

5.0 

Gd 

289 

0.61 

0.49 

0.72 

3.2 

Y 

212 

0.46 

0.48 

0.77 

4.9 

Yb 

148 

0.21 

0.47 

0.77 

4.1 

Tb 

0.018 

- 

0.46 

- 

- 

around  Gd  radius  with  intermediate  values  of  t  and  Vf  (Fig.  2b  and 
Fig.  9). 

The  Arrhenius  plots  of  the  macroscopic  grain-boundary  resis¬ 
tance  L/(S  x  Rgb)  shows  a  similar  trend  to  the  bulk  one,  reaching 
a  maximum  for  Gd-doped  samples  (Figs.  9  and  10).  The  activation 
energies  of  the  grain-boundary  resistance  are  larger  than  those  cor¬ 
responding  to  the  bulk,  ranging  from  0.72  to  0.82  eV  and  reaching 
the  minimum  value  for  Gd-doped  sample  (Table  2).  The  different 
grain  boundary  resistances  cannot  be  only  attributed  to  differences 
in  grain  sizes,  since  the  average  grain  size  is  similar  for  all  the  sam¬ 
ples,  ranging  between  1  and  2  p,m.  It  should  be  considered  that  the 
grain  boundary  conductivity  calculated  using  Eq.(4)  is  not  reliable 
because  the  thickness  of  the  grain  boundary  parallel  to  the  current 
flow  is  much  smaller  than  that  of  the  pellet  (L/S).  Thus,  both  the 


grain  boundary  layer  thickness  and  grain  size  must  be  considered 
to  evaluate  the  microscopic  grain  boundary  conductivity,  defined 
as  [39,40]: 


<^gb 

D 


(6) 


where  D  is  the  average  grain  size  and  <$gb  the  thickness  of  the  grain 
boundary  layer.  The  grain  boundary  layer  widths  8gb  and  the  grain 
size  D  should  be  known  to  obtain  <7gb  and  they  can  be  only  obtained 
by  microscopic  analysis.  An  alternative  approach  to  estimate  the 
grain  boundary  conductivity  consists  in  the  use  of  the  bulk  and 
grain  boundary  resistances  and  capacitances  [40]: 


CTgb  —  Wgb£?o£gb 
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Fig.  9.  Variation  of  the  bulk  conductivity  and  the  macroscopic  graing  boundary  (GB) 
resistance  at  150°C  as  a  function  of  the  rare-earth  ionic  radius. 


The  grain  boundary  width  5gb  can  be  also  estimated  using  the 
Eq.8,  which  has  been  obtained  by  combining  Eqs.  (6)  and  (7)  and 
considering  a  similar  value  of  permittivity  for  grain  interior  and 
grain  boundary  (£gb  ^  £b): 


£gb  -  c 


£0£b  DS 


gb 


(8) 


The  estimated  grain  boundary  widths  vary  between  3.2  and 
4.5  nm,  which  are  comparable  to  the  values  reported  for  other  solid 
electrolytes  (e.g.  l-2nm  for  Ceo.sGdo^Oi.g  [41]  and  3-4 nm  for 
Zro.84Yo.i60i.92  [42]).  The  grain  boundary  widths  in  Table  2  have 
been  calculated  presuming  a  similarity  of  permittivity  between  the 
grain  boundary  and  the  grain  interior  regions,  although  the  com¬ 
position  of  the  grain  boundaries  might  be  different  depending  on 
segregations  and  morphology  changes,  so  that  the  values  of  <$gb  are 
only  estimated,  but  they  are  reasonable  when  compared  to  other 
solid  electrolytes. 

The  Arrhenius  plots  of  the  total  conductivity  in  wet  air  and  5% 
H2-Ar  gas  mixture  are  shown  in  Fig.  11.  Similarly  to  the  bulk  con¬ 
ductivity,  the  total  conductivity  depends  on  the  earth-rare  dopant 
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Fig.  10.  Arrhenius  plots  of  the  grain  boundary  resistance  for  BaCeo.9Ln0.i03_5  samples  under  (a)  wet  air  and  (b)  wet  5%  H2-Ar. 
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Fig.  11.  Arrhenius  plots  of  the  total  conductivity  for  BaCeo.gLno.iO^  samples  under  (a)  wet  air  and  (b)  wet  5%  H2-Ar. 


size,  the  Gd-doped  samples  showing  the  highest  conductivity  val¬ 
ues  (Fig.  12). 

It  should  be  commented  that  Bonanos  et  al.  [18]  reported  pre¬ 
viously  larger  conductivity  values  in  Gd-doped  BaCe03  than  in 
Y-doped  BaCe03. 

It  is  also  evident  from  Fig.  1 1  that  doped  cerates  show  higher 
conductivity  in  air  than  in  hydrogen  atmosphere  at  the  high  tem¬ 
perature  range.  In  air  with  high  oxygen  partial  pressure,  these 
specimens  show  significant  p-type  conduction  and  when  water 
vapour  or  hydrogen  is  introduced  into  the  atmosphere,  the  p-type 
contribution  decreases  compared  to  proton  contribution. 

The  conductivity  of  all  these  electrolytes  is  substantially  higher 
than  that  of  the  ceria-based  ones  (e.g.  Ceo.sGdo.2O1. 9)  with  oxide 
ion  conduction  and  especially  at  temperatures  lower  than  500  °C 
(Fig.  11a)  and  therefore  they  are  good  candidates  for  solid  oxide 
fuel  cell  applications. 

The  conductivity  of  Gd-doped  sample  as  a  function  of  the  oxy¬ 
gen  partial  pressure  is  shown  in  Fig.  13a.  The  water  partial  pressure 
was  kept  at  0.02  atm  by  bubbling  water  at  room  temperature.  At 
the  high  p02  range  of  10_5-0.21  atm,  the  conductivity  increases 
with  p02,  indicating  that  the  conduction  mechanism  includes 
ionic  and  electron-hole  components.  A  plateau  region  observed 
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Fig.  12.  Variation  of  the  total  conductivity  at  300, 500  and  600  °C  in  BaCeo.gLno.i  03_^ 
compounds  as  a  function  of  the  rare-earth  radius  of  the  dopant  in  wet  air  atmo¬ 
sphere. 
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Fig.  13.  Oxygen  partial  pressure  dependence  of  the  total  conductivity  at  700  and 
800  °C  for  (a)  BaCeo.gGdo.iO^  and  (b)  BaCe0.9Tb0.i  03_5  samples.  Solid  lines  corre¬ 
spond  to  the  fitting  results  using  Eq.  (9). 

at  p02  values  lower  than  10-5  atm  is  attributed  to  ionic  conduc¬ 
tion.  Consequently,  the  total  conductivity  in  these  conditions  can 
be  approximated  by  the  classical  model: 

at  =  o-,+CTp(p02)1/4  (9) 

where  crz  and  cr°  are  the  values  of  ionic  and  partial  p-type  electronic 
conductivities  at  unit  oxygen  pressure.  This  model  was  used  for 
fitting  the  o  vs.  p02  curves  (Fig.  13a). 

Tb-doped  sample  exhibits  a  different  dependence  of  conductiv¬ 
ity  with  p02  (Fig.  13b).  In  this  case,  the  total  conductivity  increases 
nearly  linear  with  decreasing  p02  with  a  slope  of  about  -1/5. 
This  behaviour  is  attributed  to  the  variation  in  oxygen  vacancy 
concentration  with  p02  due  to  gradual  reduction  of  Tb4+  to  Tb3+. 
Consequently  the  ionic  conduction  is  not  constant  and  the  total 
conductivity  does  not  follow  the  simple  model  of  (Eq.  (9)). 

4.  Conclusions 

A  freeze-drying  method  was  used  to  prepare  barium  cerates 
powders  with  different  dopants  on  the  B-site  of  the  perovskite 
structure  BaCe0.9Ln0.iO3_5  with  Ln  =  La,  Nd,  Sm,  Gd,  Yb,  Tb  and 
Y.  All  materials  crystallize  in  the  same  Pmcn  space  group.  The 
lattice  cell  volume  decreases  as  the  ionic  radius  of  the  dopant 
decreases  as  expected  for  a  solid  solution.  Only  Tb-doped  sam¬ 
ple  shows  a  different  behaviour  due  to  the  presence  of  mixed  Tb3+ 
and  Tb4+  ions  in  the  structure  as  confirmed  by  XPS  analysis  and 
impedance  spectroscopy.  Goldschmidt  tolerance  factor  and  free 
volume  indicated  that  Gd  is  the  best  dopant  for  optimizing  oxy¬ 
gen  ion  conductivity  in  BaCe03 -based  materials.  Dense  ceramic 
pellets  with  similar  microstructure  for  all  materials  were  obtained 


after  sintering  at  1400  °C.  The  main  difference  in  the  microstruc¬ 
ture  is  the  average  grain  size,  which  decreases  as  the  ionic  radius 
of  the  dopant  is  getting  smaller.  The  bulk,  grain  boundary  and  total 
conductivities  of  these  materials  were  studied  by  impedance  spec¬ 
troscopy  as  a  function  of  the  rare-earth  doping.  The  values  of  the 
bulk  and  grain  boundary  conductivities  vary  with  the  rare-earth 
size,  reaching  a  maximum  for  Gd-doped  samples.  This  result  con¬ 
tradicts  the  commonly  accepted  belief  that  Y-doping  is  the  best 
option  for  optimising  the  ionic  conduction  in  BaCe0.9Ln0.iO3_5  sys¬ 
tem. 
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